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ABSTRACT 

We examine the nature of RXJ0806. 3+1527 and show that it is possible to reconcile the observed period 
decrease and X-ray luminosity with the transfer of mass between two white dwarfs provided that: either the 
system is (i) still in the early and short-lived (< 100 yr) stages of mass transfer due to atmospheric Roche-lobe 
overflow, or (ii) in a standard, long-term, quasi-stationary mass-transfer phase that is significantly (~ 90%) 
non-conservative and the conversion of accretion energy to X-rays is quite inefficient. In either of the two cases 
and for a wide range of physical parameters, we find that orbital angular momentum is lost from the system at 
a rate that is a factor of a few (< 4) higher than the rate associated with the emission of gravitational waves. 
Although the physical origin of this extra angular momentum loss is not clear at present, it should be taken into 
account in the consideration of RX J0806. 3+1527 as a verification Galactic source for LISA. 
Subject headings: Stars: Binaries: Close, Stars: White Dwarfs, Gravitational Waves 



1. INTRODUCTION 

In recent years, much attention has been devoted to the na- 
ture of the X-ray source RXJ0806. 3+1527. The source was 
discovered by ROSAT in 1990 (Beuermann et al. 1999) and 
found to be variable with a period of ~ 321 s by Israel et al. 
(1999). The latter authors tentatively ascribed the periodic- 
ity to the rotation of a white dwarf in an intermediate polar 
type cataclysmic variable (see also Beuermann et al. 1999; 
Burwitz & Reinsch 2001; Norton, Haswell, & Wynn 2004). 
The absence of a second period (corresponding to the orbital 
motion) and the presence of He emission lines in an optical 
follow-up, however, lead Israel et al. (2002) to reclassify the 
system as an AM C Vn type double degenerate. In this picture, 
the ~ 321s periodicity reflects the orbital motion, making 
RXJ0806. 3+1527 the tightest known binary to date (see also 
Burwitz & Reinsch 2001; Ramsay, Hakala, & Cropper 2002) 
and an excellent verification source for the gravitational-wave 
space mission LISA planned by NASA and ESA. 

More recent observations revealed two possibly major 
problems for the AM CVn model: (i) Israel et al. (2003), as- 
suming a distance of 500 pc, derived a low average X-ray lu- 
minosity of ~5 x 10 32 ergs _1 , and(ii) Hakala etal. (2003)and 
Strohmayer (2003) found the observed period to be decreas- 
ing with time at a rate of ~ 10~ n ss" 1 (see also Israel et al. 
2004; Strohmayer 2005). The X-ray luminosity poses a prob- 
lem because it is significantly below the luminosity expected 
from mass accretion in compact double white dwarfs, while 
the period decrease poses a problem because mass transfer in 
AM CVns is expected to expand the orbit rather than shrink 
it. A third possible model was therefore proposed to be the 
unipolar inductor model developed by Wu etal. (2002). In this 
model the X-ray flux is generated by electric currents between 
two close but detached white dwarfs. However, shortcomings 
of this model have recently been pointed out by Barros et al. 
(2005) and Marsh & Nelemans (2005). The question on the 
nature of RX J0806. 3+1527 therefore still remains open. 

In this Letter, we revisit the AMCVn model for 
RXJ0806. 3+1527 and investigate under which conditions the 
assumption of mass transfer (MT) between two white dwarfs 
can comfortably explain all the currently available observa- 
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tional constraints. Our motivation stems from both trying to 
understand the origin of the X-ray emission and the measured 
period decrease and to examine the physical properties of this 
binary and its role as a verification source for LISA. 

2. THE MASS AND RADIUS OF THE DONOR STAR 

We assume RXJ0806. 3+1527 to be an AMCVn type dou- 
ble degenerate with an orbital period P or b = 5.4 min (~ 321 s). 
We denote the masses of the two white dwarfs (WD) by Mi 
and Mi, and assume M\ > M2 so that the primary corresponds 
to the mass accretor and the secondary to the mass donor. 
Since the evolutionary channels leading to the formation of 
AM CVns involve multiple MT phases (e.g., Han 1998, Nele- 
mans et al. 2001a), the orbit can be safely assumed circular. 

The assumptions that the observed periodicity of 5.4 min 
corresponds to the orbital period and that the secondary fills 
its critical Roche lobe can be used to constrain the mass M2 
and radius R 2 of the donor. For this purpose, we first deter- 
mine the radius of the secondary as a function of its mass by 
means of the mass-radius relation for zero-temperature WDs 
(Nauenberg 1972): 
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where M c h is the Chandrasekhar mass (see also Han & Web- 
bink 1999). The resulting radii are shown in Fig.^solid line). 
Next, we consider the mass-radius relations for finite temper- 
ature helium WDs with masses from 0. 1 to 0.4 M Q derived by 
Hansen & Phinney (1998). Radii for effective temperatures 
T 2 = 5 x 10 3 and 10 4 K, and hydrogen envelopes of M H = 10" 6 
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and3x 10 _4 M Q are shown in Fig.[T](long-dashed lines). They 
are systematically larger than radii for zero-temperature WDs, 
with the largest differences occurring at the lowest masses. 

Under the assumption that the secondary is not signifi- 
cantly out of thermal equilibrium, its radius is approximately 
equal to the volume-equivalent radius R^ 2 of its critical Roche 
lobe (see Eggleton 1983). The variations of the latter as a 
function of M 2 are shown in Fig. ^ (short-dashed lines) for 
P otb = 5.4 min and Mi = 0.2,0.6, 1.0, 1.4 M . It follows that, 
for zero-temperature WD models, R 2 = R^ 2 forM 2 = 0.1 3M Q 
and R 2 = 0.024 ,r? Q (cf. Israel et al. 2002). Finite-temperature 
mass-radius models yield larger donor masses and radii. For 
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FIG. 1. — Mass-radius relations for zero-temperature (solid line) and finite 
temperature (long-dashed lines) WD models. From left to right, the finite 
temperature models correspond to (Mu/Mq^/K) = (10 -6 ,5 X 10 3 ), (3 X 
10^,5 x 10 3 ), (KT 6 ,10 4 ), (3 x KT 4 ,10 4 ). The donor Roche-lobe radii 
for P orb = 5.4 min, M, = 0.2,0.6, 1.0, 1.4M , and M 2 < Mi are also shown 
(short-dashed lines). 
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FIG. 2. — Variations of Lx as a function of M\, for zero-temperature WDs, 
Mi = 0.13 Mq , a. = 0.1, and /3 = 0. 1 or 1 . Solid lines correspond to the quasi- 
stationary MT rate with J a! \, = 7qr [Eq. |3j]; dashed lines to the turn-on rate 
[Eq. 1411. The grey horizontal band indicates the observationally inferred Lx 
for distances of 100-1000 pc. Accretors with M[ < 0.13Mq (hatched area) 
are excluded because of the Mi > M2 assumption. 



T 2 = 10 4 K and M H = 3 x 1O _6 M , e.g., M 2 = 0.24 M Q and 
R 2 = 0.029 R Q . 

3. THE MASS-TRANSFER RATE AND X-RAY LUMINOSITY 

In AM C Vn binaries, the observed X-ray luminosity is gen- 
erated by mass accretion onto the primary. Ramsay et al. 
(2005) have shown that a large fraction of the accretion lumi- 
nosity in AM C Vns is radiated in the UV rather than in X-rays. 
In particular, they estimated the ratio of the X-ray luminosity 
to the UV luminosity to be - 1 /1000 for AM CVn and ~ 1 /3 
for GP Com. In order to accommodate this, we set Lx = aL icc , 
where Lx and L acc are the X-ray and total accretion luminos- 
ity, respectively, and < a < I. We conservatively set a = 0. 1 
and we note that a can very well be lower by 1-2 orders of 
magnitude. The linear dependence of L x on a allows for an 
easy rescaling of our results for different values of a. 

The accretion luminosity due to MT is given by 

iacc=-/3M 2 ($Ll-$Rl), (2) 

where (3 is the fraction of the transferred mass accreted by the 
primary, and $li and $ri are the Roche potential values at 
the inner Lagrangian point and at the surface of the accretor, 
respectively (see Han & Webbink 1999 for details, but note 
the reverse definition of M\ and M 2 ). To determine $m, we 
model the accretor as a zero-temperature white dwarf, noting 
that finite temperatures lead to an increase in the radius and 
thus a decrease in the absolute value of <£>ri and L acc . 

3.1. Quasi-Stationary Mass Transfer 

In order to estimate the accretion rate and X-ray luminosity, 
we need to determine M 2 from the donor. We assume that MT 
turns on instantaneously when R 2 = Ri 2 and that M 2 is given 
by the quasi-stationary rate 

A-rVfir-'R 0) 

t,s2-(,L2 \"2 -/orb / 

Here 7 or b is the orbital angular momentum, R 2 and / orD cor- 
respond to time derivatives of R 2 and J m ±, at constant mass, 
Cs2 = (inR 2 /lnM 2 ) s is the donor's adiabatic radius-mass ex- 
ponent, and Cl2 = (lnRL 2 /lnM 2 ) (e.g. Rappaport, Joss, & 



Webbink 1982; Ritter 1988; Han & Webbink 1999). Expres- 
sions for Q 2 and £l2 appropriate for WDs are given by Han & 
Webbink (1999) 1 . For our purpose, we neglect changes in the 
stellar radius related to the cooling of the white dwarf at con- 
stant mass, and assume the orbital evolution in the absence of 
MT to be driven entirely by orbital angular momentum losses 
due to gravitational radiation (GR). In Eq. yj, we then have 
R 2 = and / orb = J CR . 

For zero-temperature WDs (i.e. M 2 = 0.13M Q ), we calcu- 
late M 2 — 1O" 7 -1O" 5 M yr" 1 for conservative MT (f3 = 1), and 
lO^-lO^MQyr -1 for non-conservative MT {(3 = 0.1). The 
associated X-ray luminosities are shown in Fig.|2]as functions 
of Mi (solid lines). For (3= 1, Lx is well outside the grey hori- 
zontal band indicating the observational constraints on Lx for 
distances between 100 and lOOOpc (Israel et al. 2002, Marsh 
& Nelemans 2005). For /3 = 0.1, the calculated Lx is consis- 
tent with the observed one if 0.13 < Mj /M Q < 0.35, a very 
narrow range of accretor masses 2 . For a given value of 0, 
finite temperature WDs give systematically higher MT rates 
and thus higher Lx values. 

We conclude that the quasi-stationary MT rate [Eq. 0] 
yields Lx values that are only marginally compatible with 
the observed range, unless a ~ 1/1000 (as for AM CVn) 
and (3 ~ 0.1. However, it should be noted that the quasi- 
stationary rate is only representative of the long-term average 
MT rate. Marsh & Nelemans (2005) estimated that deviations 
from the quasi-stationary rate can take place on time scales of 
~ 100 yrs, suggesting that below-average MT rates could be 
sustained for time spans sufficiently longer than the 15 years 
since the discovery of the source. The mechanism responsi- 

1 The £ s 2 expression derived by Han & Webbink (1999) is valid only 
for zero-temperature WDs [M2-R2 relation as in Eq. GJ]. In our finite- 
temperature WD calculations we therefore use a generalization of f s 2 based 
on the M2-R2 relations of Hansen & Phinney (1998). The adopted expression 
for £l2 is derived under the assumption that any mass lost from the system 
carries away the specific orbital angular momentum of the accretor. 

2 Equation |2j is often approximated by L acc = -/3GMiM 2 /Ri- However, 
as noted by Han & Webbink (1999), this is inappropriate for compact double 
WDs in which the donor is located deep within the potential well of the accre- 
tor. For the system parameters considered here, the use of the approximation 
yields Lx values up to a factor of ~ 5 larger than those shown in Fig. [2] 
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ble for driving the MT rate away from its equilibrium value 
remains unclear, but, as noted by Marsh & Nelemans (2005), 
a comforting list of candidates exist, although their applica- 
bility to double WDs has yet to be investigated. 

3.2. Atmospheric Roche-lobe Overflow 

Alternatively, the MT rate in RX J0806.3+1527 could also 
be below its long-term average if the system started MT only 
recently, and therefore is in a stage of atmospheric Roche- 
lobe overflow (RLO) and has yet to evolve towards the quasi- 
stationary state. A modified rate applicable to turn-on MT 
was derived by Ritter (1988) who pointed out the role of the 
atmospheric properties of the donor star during the turn-on 
phases (see also D' Antona, Mazzitelli, & Ritter 1989; Kolb & 
Ritter 1990). We therefore reconsider the Lx derivation under 
the assumption that MT started only recently. 

Following Ritter (1988), we determine 
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with an associated exponential turn-on time scale 

Ml 2 R 2 / orb 

Here G is the Newtonian constant of gravitation, 7Z the uni- 
versal gas constant, T 2 the effective temperature, p 2 the pho- 
tospheric mass density, [i 2 the photospheric mean molecular 
weight, and Hp 2 the photospheric pressure scale height. Fur- 
thermore, F 2 (q) is a function of the mass ratio q = M\/M 2 
which, for 0.5 < q < 10, can be approximated as F 2 (q) = 
1.23 + 0.5 log q. 

From Eq. it is clear that, although convenient for the 
derivation of mass-radius relations, the zero-temperature ap- 
proximation is not suitable when considering the turn-on MT 
rate in semi-detached double WDs. Instead we determine M 2 
for a WD donor with a finite but low effective temperature 
T 2 = 2500 K. We set the mean molecular weight p 2 = 4, as 
appropriate for a non-ionized helium-rich gas, and adopt the 
equation of state of an ideal gas to determine the mass den- 
sity p 2 from the pressure P 2 and temperature T 2 . From the 
temperature-pressure stratifications for cool white dwarfs de- 
rived by Saumon & Jacobson (1995) and Rohrmann (2001), it 
then follows that P 2 ~ 10 9 dyncm~ 2 and p 2 ~2x 10~ 2 gcm~ 3 . 
The corresponding pressure scale height for M 2 = 0.13M Q 
and R 2 = 0.024 R Q is ff K ~5x 10~ 6 R 2 . Setting R 2 = R 
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then yields MT rates M 2 ~ 5 x 10 M yr and exponential 
turn-on time scales ~ 2-12 yrs. The corresponding X-ray 
luminosities are shown in Fig.[2](dashed lines). They are con- 
siderably lower than those obtained using the quasi-stationary 
MT rate and are consistent with the observed Lx for a wide 
range of M\ values, especially for (3 ~ 0.1. We note however 
that if a ~ 1 / 1000, the j3 ~ 1 case is more favorable. 

The derivation for WDs of 5 x 10 3 K and 10 4 K yields 
M 2 ~ 2 x 10" 10 and ~ 10" I2 M© yr" 1 , and turn-on time scales 
of ~ 5-15 and ~ 5-20 yrs, respectively. The associated Lx 
are therefore even lower than those shown in Fig.|2] 

4. THE ORBITAL ANGULAR MOMENTUM LOSS RATE 

The AM CVn model for RX J0806.3+1527 implies that the 
orbital period should be increasing in time, unless additional 
orbital angular momentum loss mechanisms besides GR are 
operating in the system; or if MT started so recently that 
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FIG. 3. — The orbital angular momentum loss rate due to spin-orbit cou- 
pling required for the AM CVn model to be compatible with the observed 
period decrease in RXJ0806. 3+1527. Solid lines correspond to the quasi- 
stationary MT rate for zero-temperature WDs [Eq. dashed lines to the 
turn-on rate for T2 = 2500 K WDs [Eq. Ell. In the latter case, the curves as- 
sociated with {3 = 0. 1 and 1 are almost indistinguishable because of the small 
contribution of the M2 terms in Eqs. 151 151 Both the zero-temperature and 
T2 = 2500 K calculations are based on a donor mass of 0. 1 3 Mq . The hatched 
region is excluded for the same reason as in Fig.l2l 



the accompanying orbital expansion has yet to overtake the 
orbital contraction driven by GR. Exploration of the latter 
possibility requires detailed MT calculations using realistic 
WD models, which is beyond the scope of this investigation. 
Instead, we focus on the additional orbital angular momen- 
tum losses required to reconcile the predicted orbital evolu- 
tion with the measured period decrease of 3.6 x 10 _11 ss _1 
(Strohmayer 2005). We tentatively ascribe the additional 
orbital angular momentum losses to spin-orbital coupling 
through tides and/or magnetic fields, although our analysis is 
generally valid for other mechanisms as well. 

From the definition of J OI \, and Kepler's third law, it follows 
that 
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Writing 7 or b as the sum of the contributions from gravitational 
radiation (GR), spin-orbit coupling (SO), and systemic mass 
loss due to non-conservative mass transfer (MT), on the other 
hand, yields 

■^orb =JgR + JsO+JmJ- (7) 

The first two terms in Eq. Q operate in both detached and 
semi-detached double WDs. The third term becomes active 
only after the onset of MT. We model this term assuming that 
any systemic mass loss carries away the specific orbital angu- 
lar momentum of the accretor, i.e. 



MT : 



(1 



orb 



-M 2 . 



(8) 



Mi M\+M 2 

Substitution of Eqs. (|6j and (|8} into Eq. @, and use of 
the well known expression for 7gr (e g Landau & Lifshitz 
1962) and the observed period and rate of period decrease 
yields the /so values required for consistency with obser- 
vations. For zero-temperature WDs (M 2 = O.DM0), /so is 
shown in Fig. [5] in units of /gr. In the case of the quasi- 
stationary MTrate [Eq. Q], / so ~ 1.5/gr for Mi ~ 1.4M , 
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and /so > 3.5/gr for M\ < 0.35M© (the Mi range where the 
calculated Lx for a = 0.1 and = 0.1 is consistent with the 
observed one). For 5 x 10 3 K and 10 4 K WDs, /so is even 
larger. On the other hand, in the case of the turn-on MT rate 
[Eq. 0], /so < 4.5 /gr over the entire M\ range, and /so is 
negligible (~ 10" 3 /gr) for Mi ~ 1.4M . This also applies to 
finite WD temperatures because of the low MT rates during 
the onset of RLO. 

5. DISCUSSION AND CONCLUSIONS 

We have revisited the AMCVn model for 
RXJ0806. 3+1527 and confronted its theoretical predic- 
tions with the observational constraints. For an orbital period 
of 5.4 min and a donor temperature of less than 10 4 K, RLO 
from the secondary occurs for M2 — 0. 13-0. 24 M© and 
R 2 ~ 0.024-0.029/?©. Under the assumptions that 10% 
of the transferred mass is accreted by the companion, and 
that 10% of accretion luminosity is radiated in X-rays, the 
X-ray luminosity due to accretion is found to be comfortably 
consistent with the observed X-ray luminosity within the 
distance uncertainties (d < 1 kpc). 

We have examined both long-term quasi-stationary and 
short-term atmospheric MT with exponential growth. We 
found that the predicted Lx and / or b are more comfortably rec- 
onciled with the observations in the latter of the two cases. 
The lower MT rates for atmospheric RLO yield lower ac- 
cretion luminosities, which significantly alleviates the con- 
straints on the fraction of the mass accreted by the companion 
and the fraction of the accretion luminosity emitted in X-rays. 
The predicted orbital evolution is compatible with the mea- 
sured period decrease, if orbital angular momentum losses 
other than just GR (e.g. due to spin-orbit coupling) generally 
play a non-negligible role. We find that, in the case of quasi- 
stationary MT, an extra loss rate of at least a factor of ~ 1.5 
times the GR rate is required. On the other hand, in the case 
of atmospheric MT, this factor is smaller than 4.5 over the en- 
tire Mi range and vanishes altogether for Mi > 1.2M Q . The 
drawback of the atmospheric RLO hypothesis is that MT must 
have started only very recently (within ~ 5 x , i.e., at most 
~100yr). 

The standard hypothesis of quasi-stationary MT, however, 
can also become consistent with observations if MT is signifi- 
cantly non-conservative and the fraction of the accretion lumi- 
nosity emitted in X-rays (instead of UV) is as low (~ 1 /1000) 
as estimated by Ramsay et al. (2005) for AMCVn. The draw- 
back of this hypothesis is that the required additional orbital 
angular momentum losses are about 1.5-4 times the GR loss 
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rate (forM; > 0.3 M Q ). Without a detailed study of spin-orbit 
coupling through tides and/or magnetic fields it is not obvious 
whether such high loss rates are sustainable. 

Besides spin-orbit coupling, the rate of orbital angular mo- 
mentum loss can also be increased if matter leaving the sys- 
tem is assumed to carry away the specific orbital angular mo- 
mentum of the donor rather than that of the accretor. In addi- 
tion, if RXJ0806. 3+1527 is a direct impact accretor, the ab- 
sence of an accretion disk provides yet another orbital angu- 
lar momentum sink since no angular momentum from the ac- 
creted material can be directly fed back into the orbit (Marsh, 
Nelemans, & Steeghs 2004). We also note that in the case 
of quasi-stationary MT, these extra orbital angular momen- 
tum losses tend to increase M2 and Lx [see Eqs. (|2j-(|3j]. For 
Mi > 0.3 M Q the increase is smaller though than the orders 
of magnitude uncertainty in the fraction of the accretion lu- 
minosity emitted in X-rays. 

We conclude that, in view of the uncertainties in the WD 
temperatures, the MT mechanism, the fraction of the accre- 
tion luminosity radiated in X-rays, and the acting orbital an- 
gular momentum loss mechanisms, the observed X-ray lumi- 
nosity and rate of period decrease cannot be used to exclude 
that RXJ0806. 3+1527 is an AMCVn type double degener- 
ate. The standard association of accretion and X-ray emission 
does not at all appear inconsistent with the current observa- 
tions and hypotheses about the origin of the X-rays being un- 
related to accretion are not necessary. 

The present study reveals the possibility of additional angu- 
lar momentum loss rates due to spin-orbit coupling in double 
WDs and clearly indicates that the use of these interacting 
binaries in LISA detection and data analysis may not be as 
straightforward as previously thought. 
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